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Abstract 

The rational and hyperbolic Ruijsenaars-Schneider models and their non-relativistic 
limits are maximally superintegrable since they admit action variables with globally well- 
defined canonical conjugates. In the case of the rational Ruijsenaars-Schneider model 
we present an alternative proof of the superintegrability by explicitly exhibiting extra 
conserved quantities relying on a generalization of the construction of Wojciechowski for 
the rational Calogero model. 



1 



1 Introduction 



Let us consider a Hamiltonian system (M,u,h), where (M,u) is a symplectic manifold of 
dimension 2n and h is the Hamiltonian. The system is Liouville integrable if there exist n 
independent functions hi G C°°(M) (i = l,...,n) that are in involution with respect to the 
Poisson bracket and the Hamiltonian h is equal to one of the hi. A Liouville integrable system 
is called maximally superintegrable if it admits (n — 1) additional constants of motion, say 
fj G C°°(M), that are time- independent, globally smooth, and the (2n — 1) functions 

hi, ■ ■ ■ ,h n , fi, . . . , f n -i (1) 

are independent, i.e., their differentials are linearly independent on a dense submanifold of M. 
As a general reference on superintegrability, we refer to pQ. Maximal superintegrability and 
the compactness of the joint level surfaces 

h i = h° i (Vi = l,...,n) and ft = /° (Vj = 1, . . . , n - 1) (2) 

for generic constants h® and /? implies the periodicity of the generic flows of the system 
(M,u,h). Maximally superintegrable systems with periodic flows are very rare, the classi- 
cal examples being the isotropic harmonic oscillator and the negative energy sector of the 
Kepler problem. Other examples are provided by magnetic analogues and higher dimensional 
generalizations of the Kepler problem [2] . See also [31 0] for interesting classification results in 
low dimensions. 

There exists a large class of Liouville integrable systems that are maximally superintegrable 
in a rather obvious manner. These are the systems that admit action variables with globally 
well-defined canonical conjugates. Typical examples are scattering systems having configuration 
space trajectories q(t) = (qi(t), . . . , q n {f)) with large time asymptotes of the form 

Qi(t) ~ it + t v t (Vz = 1, . . . , n) as t — > oo (3) 

in such a way that the action variables pf = m^vf and their canonical conjugates belong to 
C°°(M) and together they parametrize the phase space M. This is the expected behaviour 
in many-body models of particles interacting via repulsive pair potentials. For example, the 
scattering characteristics of the rational and hyperbolic Calogero-Sutherland models El [7] 
as well as of their Ruijsenaars- Schneider (RS) deformations [S] have been analyzed in [9], and 
the results imply that these many-body models are maximally superintegrable. 

In our opinion, it is interesting to find explicit expressions for the constants of motion even 
for those systems whose maximal superintegrability is already known from abstract arguments. 
This is especially true if simple expressions can be obtained. A nice example is the direct proof 
of the maximal superintegrability of the rational Calogero model presented by Wojciechowski 
[TO] . His arguments were later generalized to the rational Calogero models based on arbitrary 
finite Coxeter groups [TT]. An attempt to study the hyperbolic Sutherland model in a similar 
manner was made in [T2], but in this case the constants of motion have a more complicated 
structure and to our knowledge fully explicit globally smooth expressions of them are still not 
available. 

The main goal of this letter is to explicitly exhibit constants of motion that show the maximal 
superintegrability of the rational Ruijsenaars- Schneider model. We find that the method of 
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Wojciechowski can be successfully applied in this case. Before explaining this in Section 3, in 
the next section we present a general argument implying the maximal superintegrability of the 
Liouville integrable systems that admit action variables with globally well-defined canonical 
conjugates. The precise definition of the foregoing condition is given below. 



2 A class of maximally superintegrable systems 

Consider the Liouville integrable systems (M,u,hi) associated with the Poisson commuting, 
independent Hamiltonians hi, . . . ,h n . Let us assume that globally well-defined action variables 
with globally well-defined canonical conjugates exist. By definition, this means that we have a 
phase space {M., fl) of the form 

M := R n x V n = {(Q, P) | Q E W 1 , P E V n } (4) 

with a connected open domain T> n C W 1 and canonical symplectic form 

n 

= ^2dPi A dQi, (5) 
i=i 

which is symplectomorphic to (M,u) and permits identification of the Hamiltonians {hi} as 
functions of the action variables {Pj}. More precisely, we assume the existence of a symplecto- 
morphism 

A : M ->■ M (6) 
such that the functions Hi := hi o A~ l do not depend on the variables {Qj} and 

yields an invertible matrix X(P) at every P e T> n . The map A is referred to as a global 
action-angle map of maximally non-compact type. 

If a global action-angle map of the above type exists, then one can introduce the functions 
fi E C°°(M) (i = 1, . . . , n) by the definition 

n n 

(f i oA- 1 )(Q,P):=Y,QjX(P)x with ^I(P) y -I(PS^,,, (8) 

3=1 3=1 

By using that A is a symplectomorphism, one obtains the Poisson brackets 

{fi, hj} M = Si j, {f, fj} M = 0. (9) 

Together with {hi, = 0, these imply that the 2n functions hi, ... , h n , fi, . . . , f n are func- 
tionally independent at every point of M. The choice of any of these 2n functions as the 
Hamiltonian gives rise to a maximally superintegrable system. For example, the (2n— 1) inde- 
pendent functions hi, ... , h n , f 2 , . . . , f n (resp. hi, ■ ■ ■ , h n , fi, . . . , f n ) Poisson commute with hi 
(resp. with /i). 
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Let us now study a Hamiltonian H G C°°(M) of the form H = Sj(hx, . . . , h n ) with some 
$3 G C°°(R n ). Under mild conditions, this Hamiltonian is also maximally superintegrable. For 
example, suppose that there exists an index I G {l,...,n} such that diSj ^ generically. 
For I = 1, the assumption guarantees that H,I%2, ■ ■ ■ , h n are independent, i.e., H is Liouville 
integrable in the sense specified at the beginning. Next, let us choose (n — 1) smooth maps 
V a ■ R n — > R n (a = 1, . . . , n — 1) that satisfy the equations 

■n 

J2v a k d k * = o (10) 

k=l 

identically on R n , and their values yield independent R n - vectors generically. In other words, the 
vectors V a (x) span the orthogonal complement of the gradient of F) at generic points igM" 
Such maps V a always exist, since the independence of their values is required only generically. 
By using V a , we define the function F a G C°°(M) by 

n 

F a := fkVaih, ...,h n ), Va = 1, . . . , n - 1. (11) 

k=l 

It is easily seen that the globally smooth functions H, h 2 , ■ ■ ■ , h n , Fi, . . . , F n _i are independent 
and they Poisson commute with H. This demonstrates that H is maximally superintegrable. 
Incidentally, the set h\, . . . , h n , Fx, . . . , F n _i also gives (2n — 1) independent constants of motion 
for H, and this holds even if we drop our technical assumption on the gradient of G C°°(R n ). 

The above arguments are quite obvious and are well known to experts (see e.g. [T3j ) . Among 
their consequences, we wish to stress the fact that all Calogero type models possessing only scat- 
tering trajectories are maximally superintegrable. Indeed, for these models action-angle maps of 
maximally non-compact type were constructed by Ruijsenaars [9]. Concretely, this is valid for 
the rational Calogero model, the hyperbolic Sutherland model, and for the relativistic defor- 
mations of these models due to Ruijsenaars and Schneider. The construction of the pertinent 
action-angle maps relies on algebraic procedures, but fully explicit formulae are not available. 
Therefore it might be interesting to display the maximal superintegrability of these models by 
alternative direct constructions of the required constants of motion. 

3 Constants of motion in the rational RS model 

Let us recall that the phase space of the rational Ruijsenaars-Schneider model is 

M = C n xR n = {(q,p)\qeC n ,peR n } (12) 

where 

C n := {q G R n | qx > q 2 > • • ■ > <?„}• (13) 
The symplectic structure u = Ylk=i ^Pk A dqk corresponds to the fundamental Poisson brackets 

{qi,Pj}M = {qu Qj}M = {Pi,Pj}M = 0. (14) 
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The commuting Hamiltonians of the model are generated by the (Hermitian, positive definite) 
Lax matrix [S] given by 



L{q,p)j,k = Uj(q,p) 
with the R + -valued functions 

Uj(q,p) := e P] JJ 



ix + (Qj ~ Ik) 



u k (q,p) 



(15) 



1 + 



r 



(16) 



where x 7^ is an arbitrary real coupling constant and the 'velocity of light' is set to unity. 
By using the diagonal matrix q := diag(gi, . . . , q n ), we define the functions Ik, I\ G C°°(M) 

by 

I k (q,p):=tr(L(q,p) k ), ll( q , p) := tr (qL(q, p) k ) V/c G Z. (17) 

It is well-known that the functions Ik pairwise Poisson commute and a convenient generating 
set of the spectral invariants of L is provided by the independent functions Ii, . . . ,I n . Note that 
one may also use the coefficients of the characteristic polynomial of L(q,p) as an alternative 
generating set. The 'principal Hamiltonian' of the model is 



h 



:(h + I-i) 



cosh(p fe ) Yl 



k=l 



1 + 



;is) 



Our subsequent considerations are based on the following important formula: 

{ J ki J j}M = jlj+k Vj, k G Z. 



(19) 



This generalizes an analogous formula found by Wojciechowski [10J in the rational Calogero 
model. The quantities Ij and I\ form an infinite dimensional Lie algebra under the Poisson 
bracket, with the l\ realizing the centerless Virasoro algebra: 



{lll)} M = {j-k)ll 



Vj, k G Z. 



(20) 



We postpone the proof of the above relations for a little while. 

Since dim(M) = In is finite, only finitely many of the functions Ik, 1^ can be independent. 
A set of 2n independent functions is given, for example, by 



T T T l T 1 

• • • j 1 n, 1 \, • • ■ , 1 n - 



To see that the Jacobian determinant 



J := det 



d(h, ■ ■ ■ , I n , I{, ■ ■ ■ , I„) 

d(p 1 ,...,p n ,qi,...,q n ) 



(21) 



(22) 



is non-vanishing generically, notice that J is the ratio of two polynomials in the 2n-variables 
e Pi ,qi (i = l,...,n), and hence it either vanishes identically or is non- vanishing on a dense 
submanifold of the phase space. It is easy to confirm the non-vanishing of J in the asymptotic 
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region where the coordinate-differences are large. In that region L(q,p) (|T5|) becomes diagonal, 
which implies the leading behaviour 



4~$> 2 ^, I^^e 3 *" (23) 



i=i 



whereby the leading term of J can be calculated in terms of Vandermonde determinants. We 
conclude from the inverse function theorem that the 2n functions f l2T|) can be serve as indepen- 
dent coordinates locally, around generic points of the phase space. 

To elaborate the consequences of ( IT9|) . let us first consider an arbitrary smooth function 
/ = . . . , I n ). Then we obtain from (TT9~j) that {{I\,I} M ,I} M = 0. This entails that I\ 
develops linearly along the Hamiltonian flow of J, 

Il(q(t),p(t)) = ll(q(0),p(0)) + t{lll} M (q(0),p(0)), Vk G Z. (24) 

In particular, the independent functions given in ( 12~T|) provide an algebraic linearization of the 
flow of / in the sense of [IT] , i.e., they can be used as explicitly given alternatives to action-angle 
type variables. 

Next, we observe that the functions 

I} M - I}m, Vj, k G Z, (25) 

Poisson commute with /. By using these, it is possible to exhibit (2n — 1) smooth functions on 
M that Poisson commute with / = . . . , /„) and are functionally independent generically. 
For instance, for any fixed j G {1,2, ... ,n} consider the functions 

C kJ := Ill 2j - Ijl k+j , ke{l,2,...,n}\ {j}. (26) 

By taking the 2n functions ( 12~T|) as coordinates around generic points of M, we can easily 
compute the Jacobian determinant 

Jj := det °$ with a,ae{l,...,n},b,Pe{l,...,n}\{j} 7 (27) 

and obtain that Jj = (l2j) n ~ 1 , which is generically non-zero. The fact that the (2n— 1) functions 
furnished by I m (m = 1, . . . ,n) and C k j ff2l)\i are independent and Poisson commute with Ij 
shows directly that Ij is maximally superintegrable. 

As another example, take the standard Ruijsenaars-Schneider Hamiltonian h in (lisp . The 
maximal superintegrability of h is ensured by the 'extra constants of motion' 

Kj := lj(I 2 - n) - - /,_!), j = 2, . . . , n. (28) 

Indeed, one has {Kj, h} M = and 

det d }n K ri = ( /2 - n ) n_1 ' a,ae{l,...,n},b,pe{2,..., n}, (29) 
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which guarantees that the (2n — 1) functions I a , Kf, are independent. Similarly, the momentum 
V := — I-i) admits the functionally independent extra constants of motion given by 
Lj := + n)- + Ij-x) for j = 2, ... , n. 

For a general / = . . . ,I n ), let us consider functions T a G C°°(M) of the form 

n 

F a :=Y,llU k a {h,...Jn), (30) 

k=l 

with some smooth maps U a : R n — > M. n (a = 1, . . . , n — 1) subject to the identity 

n n 

£(£rt«§K = o, (3D 

k=l j=l 3 

which guarantees that {J^-Oa/ = (cf. equations ffTUj) and ffTTl) ). One can always choose the 
maps ?7 a (a = 1, . . . , n — 1) in such a way that their values yield linearly independent W 1 - vectors 
at generic arguments. Then the (2n — 1) functions Ji, J-i, ... , J-^-i are independent, 

globally smooth and Poisson commute with /. 

To summarize, we have seen that the relation ffT9|) leads to an explicit linearization of the 
Hamiltonian flow associated with any / = . . . , I n ) and allows us to display the maximal 
superintegrability of the Ruijsenaars-Schneider Hamiltonian h (fTSl) in an explicit manner (and 
similarly for V and Ij for j = 1, . . . , n). The above arguments also imply, among others, the 
maximal superintegrability of any polynomial Hamiltonian / = . . . , 

Now we prove the relations (TT^|) and (TSO"]) . Direct verification is possible in principle, but it 
would require non-trivial calculations. However, it is quite easy to obtain the claimed relations 
by utilizing the derivation of the rational RS model in the symplectic reduction framework 
presented in [T3]. To explain this, we start by recalling from p3] the relevant reduction of the 
phase space 

T*GL(n, C) x 0( X ) = GL(n, C) x gl{n, C) x 0( X ) = {(<?, J R , £)}• (32) 

Here 0(x) is a minimal coadjoint orbit of the group U(n), which as a set is given by 

0(x) ■= Ox(ln - w f ) I v e C n , \v\ 2 = n}. (33) 

In (|32|) we use the trivialization of T*GL(n, C) by left-translations and identify the real Lie 
algebra gl(n, C) with its dual space with the aid of the 'scalar product' provided by the real 
part of the trace 

(X, Y) := sRtr(XF), VX, Y E gl(n, C). (34) 
In terms of evaluation functions, the not identically zero fundamental Poisson brackets read 

{g,(X,J R )} = gX, {(X,J R ),(Y,J R )} = -([X,Y],J R ), {(X + , 0, (Y + , 0} = ([X+, Y + ], 0, 

(35) 

where X + , Y + are the anti-hermitian parts of the constants X, Y G gl(n,C). The reduction is 
based on using the symmetry group K := U(ri) x U(n), where an element (i]l,i]fi) G K acts 
via the symplectomorphism ^/^ LtTlR - ) defined by 

^(vL,m) : (9, J R , ^ (Vl9Vr\ VrJ R Vr\ Vl^Vl 1 )- (36) 
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We first set the corresponding moment map to zero, in other words introduce the first class 
constraints 

J« = and ( 9 JV) + + e = 0, (37) 

and then factorize by the action of K. The resulting reduced phase space can be identified 
with the Ruijsenaars-Schneider phase space (M, uj) (IT2|) by means of a global gauge slice S. To 
describe S, we need to introduce an C(x)-valued function on M by 

€(q,p) '■= ix(ln - v(q,p)v(q, p) ] ) with v(q,p) := L(q,p)~*u(q,p), (38) 

where L(q,p) is the Lax matrix (1151) and u(q,p) is the column vector formed by the components 
Uj(q,p) ([16]) . Denoting its elements as triples according to (1321) . the gauge slice S is given by 

S := {(L(q,p)K -2q,£(g,p)) | G C n x R n }. (39) 

It has been shown in [14J that S is a global cross section of the A-orbits in the constraint- 
surface defined by (j3Tjl . and the reduced Poisson brackets (alias the Dirac brackets) reproduce 
the canonical Poisson brackets (|14[) . 

Relying on the above result, we identify (M, uj) with the model S of the reduced phase space. 
We can then realize the functions Ik,ll £ C°°(M) (TT7|) as the restrictions to 5 of respective 
A-invariant functions T k) X\ on the unreduced phase space ([52]) furnished by 

X fc (^,J R ,0:=tr((^) fe ), Z^,^0:=-^tr((^J*) VA; G Z. (40) 

Now the point is that the relation {X^,Xj} = jlj+k follows obviously from f )35|) . and this implies 
( Il9p by restriction to 5 ~ M by using that the Poisson brackets of the i^-invariant functions 
survive the reduction. The relation 

[Tlx]} = {j - k)Xl +] Vj,keZ (41) 

is also easy to confirm with the aid of (I3"5"j) . which implies (|20|) . 

We finish with a few remarks. First, we note that the variables (j2ip can be useful for 
constructing compatible Poisson structures for the rational RS model. In this respect, see [T5| 
[16] and references therein. Second, one may also construct RS versions of the Calogero constants 
of motion considered in [17] . Third, it could be interesting to study quantum mechanical 
analogues of the constants of motion and to characterize their algebras. It is very likely that 
the algebra of equations (IT9"|) . (I2U1) survives quantization. One may address this question by 
generalizing the method applied in [18] to quantize the analogous algebra in the Calogero case. 
It should be also possible to construct a quadratic algebra for the model by suitably replacing 
the Dunkl operators used in [18J with Dunkl-Cherednik operators. The pertinent algebras 
are expected to be closely related to the bispectral property of the rational RS model [T9"] . 
Finally, it is natural to ask about generalizations concerning the hyperbolic RS model and its 
non-relativistic limit. We plan to return to some of these issues elsewhere. 
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